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1 Introduction 

It is often mentioned that photons, or more generally electromagnetic probes, 
can be used to signal the formation of quark-gluon plasma (QGP) in ultra- 
relativistic heavy ion collisions. They could also be used to measure the 
temperature of this plasma at the time of its formation. Data relevant 
for this problem will soon cover a wide energy range from a nucleon-nucleon 
center-of-mass energy about 20 GeV (Pb-Pb collisions at CERN) to 200 Gev 
(Au-Au collisions at RHIC, Brookhaven) and 5.5 TeV (Pb-Pb collisions at 
LHC after 2007). 

Many production mechanisms of electromagnetic probes exist in heavy 
ion collisions. On a very qualitative level, one distinguishes the emission 
of photons in the initial stage of the collision (primary production): it is 
similar to the case of proton-proton collisions and goes via QCD Comp- 
ton scattering, quark-antiquark annihilation and bremsstrahlung emission, 
involving partons in the incoming hadrons. It is calculated in perturbative 
QCD to the next-to-leading order (NLO). The energy spectrum is power be- 
haved and it is expected to dominate the high momentum region. Hadrons 
such as TT^ and r] are also produced by perturbative mechanisms and they 
decay into photons which contribute an important background to the pho- 
ton signal at large momentum. In the heavy ion collision many partons are 
produced which eventually thermalise into a quark-gluon plasma. This hot 
QGP bubble expands and cools until a hadronic phase is formed. Photons 
are produced in the QGP phase as well as in the hot hadronic phase with 
a rate which is expected to be exponentially decreasing with the energy of 
the photon (secondary production). This thermal contribution is expected 
to be important mainly in the "intermediate" energy region region up to 10 
GeV, say. Thermally produced tt^'s and r/'s again contribute an important 
background. 

In the following we review recent progress in the calculation of photon 
rates in the QGP phase and present phenomenological applications for RHIC 
and LHC heavy ion collisions. 

2 Thermal production of hard photons 

We assume in this section that the plasma is in chemical and thermal equi- 
librium at temperature T. The rate of production of a real photon of mo- 
mentum P = {E,p), per unit time and volume, is 



where Il^^{E + ie, p) is the retarded photon polarisation tensor. The Bose- 
Einstein factor rig^E) provides an exponential damping when E ^ T. A 
similar formula holds for lepton pair production. 

The two-point correlation function is calculated in the framework of the 
hard thermal loop (HTL) resummed theory of Braaten and Pisarski 
In this approach one distinguishes two scales: the "hard" scale, typically 





of order T or larger (the energy of quarks and gluons in the plasma) and 
the "soft" scale of order gT where 5, the strong coupling, is assumed to 
be small. Collective effects in the plasma modify the physics at scale gT, 
i.e. over long distances of 0{l/gT). Propagators and vertices are modified 
by the summation of higher order diagrams. This is easily illustrated with 
the example of the fermion propagator, S{Q), which in the "bare" theory 
is simply 1/p (we neglect spin complications and make only a dimensional 
analysis)^. The thermal contribution to the one loop correction is found to 
be S(Q) ~ g'^T'^/q which is of 0{gT) when q ~ gT. The resummed propa- 
gator *S{Q) = l/{q — S(Q)) is then deeply modified, compared to the bare 
propagator, for soft momenta of 0{gT) whereas thermal corrections appear 
essentially as higher order effects for hard momenta. Thermal resummation 
affects the propagator in two ways: 1) in the space-like region [Q^ < 0), the 
propagator acquires an imaginary part, due to Landau damping (a new fea- 
ture compared to the T = theory) characterised by a thermal mass scale; 2) 
in the time-like region a thermal mass, of asymptotic value rriq = g^T^Ci?/4, 
is generated. In both cases thermal effects screen potential soft or collinear 
singularities. Likewise, the gluon propagator is modified: Debye screening 
mass = g'^T'^{Nc + Nf/2)/2> in the space-like region and quasi-particle 
masses in the time-like region.^ One-loop corrections also modify the ver- 
tices when the external momenta are soft One can construct an effective 
Lagrangian |2j in terms of effective propagators and vertices and calculate 
observables in perturbation theory. 
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Figure 1: Left diagram: one-loop contribution; right diagrams: annihila- 
tion and Compton scattering processes obtained when cutting the one-loop 
diagram. 

In the one-loop approximation, the production rate of hard photons is 
given by the diagram shown in Fig.^where the symbol • means that effective 
propagators are used. For hard photon momentum it is indeed enough to 
consider bare vertices and only one effective fermion propagator, since one 
propagator is necessarily hard and needs no resummation. The result is (for 
one quark flavor) [3] 

Imn„''(^,p) = - a g"^ T^ f ln(— ) + constant) . (2) 
9 V / 

One notices the logaritmic growth of the 2-point function with the energy of 
the photon while the thermal quark mass acts as a soft cut-off which screens 
the singularity in the forward scattering of the processes in Fig. ^ 

low case letter denotes the modulus of the 3-momentum. 
'^Nc is the number of colors, Nf the number of flavors and Cf the Casimir of the 
fundamental representation. 





Figure 2: Left diagrams: two-loop contributions; right diagrams: 
bremsstrahlung and off-shell annihilation processes obtained when cutting 
through the two-loop diagrams. 

At this point, the HTL picture of photon production, dominated by 
Compton and annihilation scattering, turned out to be at variance with the 
result of a semi-classical calculation where it was found that the production 
was dominated by the bremsstrahlung radiation of photons from quarks [1] . 
The way out of this paradox is to consider the two-loop diagrams in the 
effective theory (Fig. ^ . A dimensional analysis reveals that the dominant 
contribution comes from the kinematical region where the quark momenta in 
the loop are hard but the gluon is soft and space-like. The physical processes 
exhibited by cutting through the loop diagrams of Fig. [21 correspond to the 
bremsstrahlung emission from a quark or an anti-quark scattering in the 
plasma (as in the semi-classical approximation) and by a new process, the 
off-shell annihilation of a quark-antiquark pair where one of the (anti-)quark 
scatters in the plasma. For two flavors, the result of the calculation is very 
simple [5|: 



where the first term arises from bremsstrahlung and the second one from off- 
shell annihilation. In the general case the numerical coefficient is replaced 
by a function of the ratio of the Debye mass over the thermal quark mass [H] . 
The two-loop HTL result is not suppressed, as expected, by powers of the 
coupling g compared to the one-loop result. This is due to an enhance- 
ment factor of 0(T^/mq) ~ associated to the collinear emission of the 



One may wonder about the usefulness of the loop expansion in the HTL 
effective theory if the two-loop rate is of the same order as the one-loop 
rate. It turns out that higher loop diagrams contribute to leading order [£]. 
The physical reason is the following. Consider a virtual quark of momen- 
tum R decaying into a photon of momentum P and a quark of momentum 
Q. The photon formation time {i.e. the life-time of the virtual state) is 
Tp = r/i?^ = qr/{p{q^ + ^q)) with the transverse momentum of the quark 
measured with respect to the photon momentum. This formation time is 
long, of 0(T/mq) ~ l/g^T, and is comparable to the mean free path of 




(3) 



photon. 



the quark in the plasma given by the inverse of the quark damping rate 
7. Quark rescattering in the medium should thus be taken into account. 
This is done by resumming the self-energy diagrams on the quark propaga- 
tor going beyond the hard thermal loop approximation and including also 
the imaginary part (related to 7) generated by rescattering. In our case 
the time-like quark propagator, which in the HTL theory was + m^) 

becomes approximately +m'^ + i^qr /p): this shows that both thermal 

mass and damping rate act as competing cut-offs. Of course, to respect 




Figure 3: An example of a 
ladder diagram. 



gauge invariance, one should resum vertex corrections in the same approx- 
imation as self-energy corrections jTj, which motivates the fact that ladder 
diagrams such as the graph of Fig.^lshould be summed. Arnold, Moore and 
Yaffe showed that such diagrams indeed contribute to leading order and 
they derived a very elegant integral equation to perform the resummation. 
The imaginary part of the photon polarisation tensor takes the form: 

r+oo g2 _|_ ^2 |. 

Im n/ (P) ?a aiVcy dq[np{r) - np{q)] —^-^Re J ^^-^ q^, • f (q^, ) , 

with r = p + q and rip^q) = l/(exp(g/T) + 1), the Fermi-Dirac statistical 
weight. The dimensionless resummed quark-quark-photon vertex, f(qj,), is 
shown by the shaded vertex in Fig. ^ while q^, is proportional to the bare 
vertex. The value of f(qy) is obtained by solving the following equation 0: 



Figure 4: The integral equation resumming the ladder diagrams. 

^J{ci^) = 2ci^+g^C,Tj ^C7(l,)[f(q,+l,)-f(q,)], (4) 

where is the formation time defined above. The collision kernel is C(lj,) = 
m'^/l^{l^ + m'^) and it includes the exchange of both longitudinal and 
transverse gluons Solving for this integral equation one can obtain a fit, 
in the large E/T range 

ImU,^{P)^ag'T' M, (5) 
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Figure 5: Various contributions to the production of photons in Au-Au 
collisions at RHIC. The open dots are the predictions of decay photons 
from the non-thermal DPMJET model for heavy ion collisions, in good 
agreement with the hydrodynamic + perturbative QCD model. 



to be compared to the linear in E behaviour of Eq. Q . This is an illustration 
of the Landau-Pomeranchuck-Migdal (LPM) suppression of hard photons 
due to multiple scattering of the quark in the plasma. 

Rates for small mass lepton pairs produced at large momentum have 
also been calculated at one- loop [Si, two-loop |Hj and including the LPM 
effect (lOj. 



3 Phenomenology at RHIC and LHC 

To make predictions for photon production in heavy ion collisions, the rates 
calculated above have to be included in a hydro dynamical code which de- 
scribes the expansion and the cooling of the plasma from chosen initial state 
conditions. In the following we use the code of Ruuskanen et al. The 
rate for thermal photon production is shown by the thick dash-dotted line in 
Fig. [21 for Au-Au central collisions at RHIC and in Fig. |^ for Pb-Pb central 
collisons at LHC. Further details on the model can be found in Ref. jl3j . 
Also shown on the figures are the rates of production of decay photons from 
thermally produced resonances (mainly and 7]) (thick dashed lines) which 
dominate in the low p^, range. We also display the rates of primary photons 
produced in perturbative QCD mechanisms (solid lines) contributing in the 
large p^, range as well as the decay photons from perturbatively produced 



resonances. In these estimates account has been taken of energy loss effects 
on the jet fragmentation in the medium using the model of Ref. |14j . 
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Figure 6: Same as the above figure for Pb-Pb collisions at LHC. 
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Figure 7: Ratio of all photons over decay photons at LHC with thermal 
production (thick lines) and without thermal production. 



It is instructive to look at the ratio of "all" produced photons over decay 
photons and to compare it with the prediction for the same ratio assuming 
no thermal production. This plot quantifies the excess of direct photons 
produced by the hot hadronic medium formed in a heavy ion collision. At 
LHC, thermal effects are important but the ratio of directly produced photon 
remains only about 10%. 
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Figure 8: Thermally produced lepton pairs (dashed line) compared to per- 
turbatively produced pairs (solid line) and background from hadronic decays 
(squares) . 

A channel which should suffer less background than direct photon is 
that of small mass lepton pair at large momentum: indeed choosing a pair 
mass above the 7r° mass should eliminate a large background. The thermal 
production rates have been calculated and included in the hydrodynamical 
code • The result is shown in Fig. |HJ 



4 Further developments 

In the discussion above we have considered two independent production 
processes: emission of the photon in the initial stages of the collision and 
production in the quark-gluon plasma and hot hadronic phase. One can 
envisage "mixed" production mechanisms where a photon is emitted in the 
interaction of a hard quark or gluon (calculable in perturbative QCD) with 
a quark or the gluon of the plasma. Since the rate of hard partons is power 
behaved while thermal interactions are exponentially damped, the resulting 
rate of photon production is expected to be power behaved and to domi- 



nate over purely thermal processes at large momentum. Such examples of 
"mixed" interactions have been considered in Ref. who made simplified 
estimates of the Compton (gq — s- 7^) and annihilation {qq — s- 7^^) processes 
when one of the initial g or g is hard and the other parton is thermalised. 
The result is shown in Fig.|Ulby the dashed line to be compared with the dot- 



Figure 9: Photon pro- 
duction from "mixed" 
processes: rate from the 
Compton and annihila- 
tion processes between 
a hard parton with a 
parton in the plasma 
(dashed line) ; thermal 
photon production un- 
der the same hypothesis 
of no transverse expan- 
sion (dotted line). For 
comparison the thermal 
rate with radial expan- 
sion, used in the pre- 
vious section, is also 
shown (solid line). 



ted line which is the estimate of the purely thermal process under the same 
hydrodynamical conditions (no radial expansion). Recently B. Zakharov es- 
timated the enhanced bremsstrahlung emission of a photon by a hard quark 
traversing the plasma ^Hj. He found that, due to finite size effects, the 
usual 1/x spectrum characteristic of bremsstrahlung of a quark in vacuum 
is replaced by an approximate 1/(1 — x) spectrum for x < 1. This is illus- 
trated in Fig. ^1 which compares the fragmentation function of a quark into 
a photon in vacuum with the same in the medium taking into account only 
one rescattering in the plasma. The change of shape of the spectrum will 
lead to an enhanced production of photon by bremsstrahlung compared to 
production in the vacuum. Multiple scattering effects will somewhat reduce 
the size of the predicted enhancement pjB]. In the figures of the previous 
section we already included medium effects on bremsstrahlung production 
of photons and this lead to a suppression of photon emission compared to 
vacuum which seems to contradict the present result. In fact, there is no 
contradiction if we consider that in the latter case the photon was emitted 
by the quark outside the plasma while the former case one estimates the 
emission of the photon in the plasma. 

It would be important to have a complete and consistent estimate of the 
mixed production rates, using a realistic hydrodynamical model: in any case, 
they will increase the thermal production rates calculated in the previous 
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Figure 10: Comparison of the fragmentation function of a quark into a 
photon in vacuum (short dashed line) and in the medium (dotted line) with 
one scattering of the hard quark in the plasma. The other curves are for an 
infinite medium. 

section. 

All previous results are obtained assuming local thermal equilibrium and 
the rates are calculated per unit time. In a series of papers Boyanovsky et 
al jl7j proposed a model to explicitely account for the finite life-time of the 
plasma. They assume that photon emission starts at an initial time ti and 
ends at time tf: under this assumption, usually forbidden processes of type 
q q'y contribute since energy conservation does not hold. As a consequence 
a transient term of type alnt, where t = tf — ti, appears besides the usual 
aagt term, and it is claimed to dominate the photon production rate at LHC 
at large enough momenta. Several objections have been raised concerning 
the validity of this picture |^ ^1 EOl- In particular the hypothesis of a 
sharp switching on and off of the interactions is probably not adequate for 
the description of the physics of the problem. 
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